Substrate effects on surface magnetism of FeAV(llO) from first principles 
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Surface magnetic properties of the pseudomorphic Fe(llO) monolayer on a W(l 10) substrate are investigated 
from first principles as a function of the substrate thickness (up to 8 layers). Analyzing the magnetocrystalline 
anisotropy energies, we find stable (with respect to the number of substrate layers) in-plane easy and hard axes 
of magnetization along the [lIO]- and [001]-directions, respectively, reaching a value in good agreement with 
experiment for thick substrates. Additionally, the changes to the magnetic spin moments and the density of the 
Fe rf-states are analyzed with respect to the number of substrate layers as well as with respect to the direction 
of magnetization. With respect to the number of W(l 10) substrate layers beneath the Fe(l 10) surface, we find 
that the first four substrate layers have a large influence on the electronic and magnetic properties of the surface. 
Beyond the 4th layer, the substrate has only marginal influence on the surface properties. 

PACS numbers: 75.70.Ak, 75.70.Ri; 73.20.At 



I. INTRODUCTION 

The pseudomorphic monolayer of Fe grown on a W(llO) 
surface is very interesting from the point of view of 
magnetism. In studies of the surface magneto-crystalline 
anisotropy, the Fe monolayer on top of a W substrate has be- 
come the system of choice, since (i) the growth of the first 
Fe monolayer is pseudomorphic, (ii) the W substrate has a 
large spin-orbit coupling, and (iii) the interface anisotropy is 
the strongest ever observed. This makes the Fe monolayer on 
a W substrate a good candidate for an ab initio benchmark 
investigation of how the properties of the magneto-crystalline 
anisotropyare influenced by the substrate. 

Depending on the coverage and the surface orientation, Fe 
films on a W substrate show very different magnetic prop- 
erties. At submonolayer coverages on the (110) surface one 
observes the formation of separate islands that are nonmag- 
netic up to a coverage of 58-60%, beyond which they be- 
come ferromagnetic islands,' which may even show an out- 
of -plane magnetization until they approach the full-monolayer 
coverage.^ A single monolayer (ML) of Fe can be grown 
pseudomorphically on top of a W(l 10) surface)^ Experiments 
performed on the pseudomorphic Fe monolaye r^'^'^'^'^'^ have 
shown that^ "the prominent magnetic feature of Fe(l 10) films 
on W(llO) is a strong in-plane magnetic surface anisotropy, 
with an easy axis [110] at right angles to the bulk easy axis 
[001]." At coverages of about 1.5 ML, double-layer patches 
(sesquilayers) form with an out-of -plane magnetic easy axis^iS 
and an antiferromagnetic orde n"''^ A pseudomorphic full 
second ML has to our knowledge not been grown experimen- 
tally, but theory predicts it would also have an in-plane easy 
axis in the 110-direction.'^ Attempts of using annealing as 
a means of making the grown Fe layers full monolayers be- 
yond the first monolayer leads to a reorganization of the Fe 
in an Fe bulk-like lattice.— For continous layers of Fe, the 
easy axis of magnetization eventually changes to the in-plane 
[001] -direction it has in bulk Fc j'^''^''^-'^ the critical thickness 
of the Fe film being reported in the range between 80 and 95 A 
(95 A reported by Ref.lH, and 80-86 A reported by Ref.[ll). 

In passing, we should stress that in contrast to the ferro- 
magnetic Fe monolayer on the W(llO) surface stand results 



of recent studies (experimenta l ' ^'^'^ and ab initio) ^^i'^^'^^ of Fe 
films on the (001) surface of W, which indicate that the first 
ML of Fe is antiferromagnetic when grown on the (001) sur- 
face. Thus, the magnetic properties of the Fe monolayer is 
very sensitive to the surface orientation of the W substrate, 
and one can not make general conclusions about the magnetic 
properties of other surfaces on the basis of investigations of 
only one. 

In the present work, we study the substrate effects on the 
magnetic properties of the pseudomorphic ML of Fe(llO) 
on W(llO) from first principles. Existing ab initio calcula- 
tions performed within the full-potential linearized augmented 
plane-wave (FLAPW) method 
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the full-potential lin- 
earized muffin-tin orbital (FP-LMTO) methodr^ as well as 
pseudo-potential calculations?' confirm the above-mentioned 
experimental ferromagnetic result with an easy axis along the 
[110] -direction. At present, there is no measurable lattice- 
relaxation effect of the magnetization direction,-^ and the 
starting point for the calculations presented in this work is the 
optimized surface of Ref . |23l The main results of our work are 
magnetocrystalline anisotropy energies that are converged ab 
initio as a function of the substrate thickness, along with mag- 
netic spin- and orbital moments that consistently confirm the 
breaking of Hund's 3rd rule for the W substrate. The quality 
of the results is illustrated by explicitly showing the conver- 
gence, along with an estimate of the accuracy. 

The paper is organized as follows: in Sec.|II]we describe 
and discuss the method that we use, in particular with respect 
to the level of accuracy necessary for the calculation of sta- 
ble magnetocrystalline anisotropy energies, and in Sec.HIIlwe 
present and discuss the results of our calculations, where we 
focus the attention on the magnetic and electronic properties 
of the surface, and how the substrate influences these proper- 
ties. Finally, in Sec.|IV] we conclude. 



II. METHOD 

Since we have the flexibility that a theoretical calculation 
provides, we do not need to limit ourselves to the experimen- 
tally possible. Thus, we are able to take an approach where 



we begin with a bilayer of Fe(l 10) and W(l 10), successively 
adding W substrate layers underneath in order to investigate 
how the the sub-surface layers influences the magnetic and 
electronic properties at the surface. 

In order to facilitate the calculation of magnetocrystalline 
anisotropy energies and how they change with the num- 
ber of substrate layers, one has to perform the calculations 
with a highly accurate and mature ab initio computer code. 
The demands on accuracy imply that a full-potential code 
is necessary. Since all ab initio methods employs an "edu- 
cated guess" for the initial solution, the code must be self- 
consistently reaching a stable solution to the eigenvalue prob- 
lem of the Schrodinger equation, e.g., by use of the Rayleigh- 
Rit2s22i^ variational approach. In addition to this, the spin- 
orbit interaction must be included in the self-consistent solu- 
tion in order to produce magnetic properties that are accurate 
enough. The ab initio method of choice for our calculations is 
the full-potential (partially linearized) augmented plane-wave 
[(L)APWh-Io] with local orbitals method^ as implemented in 
the Wien2k computer code.i^ 

Within the Wien2k code, one constructs a supercell per- 
pendicular to the surface, ensuring that enough vacuum is in- 
serted between the surface and the border of the unit cell (in 
our case the vacuum amounts to at least 16 interlayer distances 
of the W(llO) substrate). In order to keep the environment 
of our surface calculation stable, we let the Fe monolayer sit 
in the middle of the unit cell of constant size, adding the W 
substrate layers under it. This choice, of course, limits our 
possibility for continously adding W substrate layers, but not 
beyond the practical limits of the code, as we shall see later 

In a study where we begin from the surface layer we expect 
to see rather big changes to the electronic and magnetic prop- 
erties as we add substrate layers to the system, due to the simi- 
larities to the properties of a quantum well. As we keep adding 
substrate layers we expect the properties to stabilize in such a 
manner that the system starts behaving like a surface layer on a 
bulk-like substrate. When this point has been reached, adding 
another substrate layer should not change the electronic and 
magnetic properties at the Fe surface layer much. 

With our demands on the accuracy, it is important that the 
surface is structurally relaxed. How a single monolayer of 
Fe(llO) relaxes on top of a W(llO) substrate has been stud- 
ied by Qian and Hiibner (Ref. |23|) . Since their results are 
very close to experiment j'-^-^'^-^'^-^''^^ and since there is no 
experimentally measurable structural effect due to changes in 
the direction of magnetization,?- we have in the present work 
chosen to adopt the optimized surface structure of Ref. |23| as 
the basis for our calculations. We have, however, in order to 
better describe the influence of the substrate on the surface 
magnetism, chosen to use a nonsymmetric layered structure 
instead of the symmetric one of Ref.|23| Thus, the structure is 
using the optimized values of the interlayer distances from the 
surface side as determined in Ref. 23 [the Fe-Wl distance is 
contracted by 12.9% and the W1-W2 distance by 0.1% com- 
pared to bulk W(l 10) interlayer distances (Wl and W2 are the 
first and second W layers under the Fe layer, respectively)], 
but has bulk-spaced W layers underneath. The Fe layer has 
an in-plane lattice misfit of about 10% compared to a pure 



Fe(l 10) surface. 

In practical ab initio calculations within density-functional 
theory, one has to make a choice of the so-called exchange- 
correlation potential, since it is not known as an exact quan- 
tity. Our calculations within the Wien2k computer code 
are performed using the [generalized gradient approximation 
(GGA)] exchan ge-c orrelation potential of Perdew, Burke, and 
Ernzerhof (Ref. l32h . In order to facilitate the calculation of 
magnetic anisotropy energies (MAEs), the variational param- 
eter, i.e., the total energy of the magnetic configuration,^^ was 
converged to an accuracy that in most cases is better than 10 
jjRy (for a detailed list of relevant accuracies, consult the re- 
sults section, in particular Tab. IJi. Also, the fluctuations in 
the charge density within the unit cell has been converged to 
less than 10^^ e/Bohr^. Since the abovementioned two con- 
vergence criteria do not automatically ensure that the eigen- 
states and eigenfunctions are sufficiently accurate, one has in 
addition to make sure that also the the calculation converges 
with respect to the parameters that control the accuracy of the 
calculations of the Wien2k code [in particular the (kinetic 
energy) cutoff value of the otherwise infinite plane-wave ba- 
sis, and the sampling in k-space]. As we shall illustrate below, 
these parameters were adjusted (increased) until the magnetic 
anisotropy energies for an increasing number of W layers sta- 
bilized. 

The parameter used to control the kinetic energy cutoff of 
the plane-wave basis depends on the muffin-tin radius used for 
the atomic part of the (L)APW basis set as follows: 

71- = ^MT^max, (1) 

where Tc is the cutoff parameter, Rmj is the muffin-tin ra- 
dius (2.35 Bohr in our calculations), and k^-^^ corresponds to 
the plane-wave cutoff (in Ry) in pseudopotential calculations, 
e.g., the value = 9 gives a k-^^^ = 199.56 eV. Thus, Tc deter- 
mines the matrix size of the eigenvalue problem, and higher 
values of Tc potentially decreases the accuracy of the resulting 
eigenfunctions and -energies. 



m. RESULTS AND DISCUSSION 

When one begins with an Fe/W(l 10) bilayer and thereafter 
add W(llO) substrate layers (on the W side of the bilayer), 
one expects that as long as the number of W layers (called Nw 
in the following) is sufficiently small, the system behaves in a 
quantum-well-like fashion, i.e., the changes to properties such 
as the magnetocrystalline anisotropy energy or the density of 
states will be large each time a substrate layer is added. At 
some critical value of Nw, one would expect that the substrate 
starts behaving like a bulk substrate and, thus, that the Fe layer 
becomes more surface-like in its properties. Below, we dis- 
cuss how the magnetocrystalline anisotropy energies and the 
electronic structure change as a function of Nw ■ 



2 



TABLE I: (Color online) Evolution of the magnetocrystalline anisotropy energy of the FeAV(l 10) surface as a function of (i) the number of 
W layers beneath the pseudomorphic Fe monolayer included in the calculation, Nw, (ii) the kinetic energy cutoff of the plane-wave basis, 
Tc = RMjkmax [below, Tc = 7 corresponds to /t^^^ = 120.72 eV, T,. = 8 to yt^^^ = 157.68 eV, and = 9 to yt^^^ = 199.56 eV], and (iii) the 
number of k-points in the full Brillouin zone, #k. The leftmost half of the table shows how the total energy of the surface magnetized along the 
easy axis, M||[110] evolve as a function of Tc, #k, and A'ly. The first column is the number of W layers, the second column is the total energy 
at {rc,#k} = {7,441}, columns 2-6 show the evolution of this total energy as the accuracy of the calculation increases (as the difference in 
total energy with respect to the previous calculation), and in column 7 is shown the total energy of the converged calculation along with the 
uncertainty estimate, as given by Eq. In columns 8-12 is shown how the MAE of the perpendicular magnetization direction (M||[110]) 
changes during convergence, where in column 12 the uncertainty estimate has been added. Column 13 lists the converged MAEs for M|| [001] 
with their respective uncertainty estimates (for the uncertainties of M 1 1 [00 1 ] relative to M 1 1 [ 1 1 0] , see the text) . Each and every number in this 
table have been determined self-consistently. 
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-3.20 


-638 
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A. Magnetocrystalline anisotropy energies 

Searching for magnetocrystalline anisotropy energies on 
the basis of ab initio methods that minimizes the total energy, 
one is looking for energy differences in the meV-range which 
has to be calculated from differences between total energies 
in the MeV-range. In order to provide us with two significant 
digits in the meV-range it is required that the computed total 
energies are numerically accurate to (within the limits of the 
physical model) the 1 1th decimal. Todays standard arithmetic 
precision (64-bit numbers with a 52-bit mantissa) provides 16 
significant digits, and thus already a sum over 1000 k-points 
reduces precision to 13 significant digits at best. Carefully 
written self-consistent full-potential methods may limit the 
further loss of precision caused by sampling, but the amount 
of k-points is limited if one wants to obtain meaningful re- 
sults. 

In order to be able to draw any conclusions about the sub- 
strate effects and the properties of the surface for a fixed 
geometry, the calculation of the MAEs must first converge 
with respect to, in particular, the kinetic energy cutoff of 
the plane-wave basis, as well as the sampling in k-spacci^ 
Since MAEs are calculated here as differences of "total en- 
ergies", this requirement automatically applies to the "total 
energies" (even though the "total energies" are not necessar- 
ily physically meaningful quantities themselves). Thus, the 
evolution of the magnetocrystalline anisotropy energies of the 
Fe/W(110) surface is described in Tab. U as a function of (i) 
the number of W layers added underneath the pseudomorphic 
Fe monolayer, (ii) the kinetic energy cutoff of the plane-wave 
basis, as defined by Eq. [T] and (iii) the number of sampling 
points (k-points) used in the full Brillouin zone (FBZ), here 
called #k in a short notation. In order to make the numbers 



directly comparable and to avoid pure k-point effects, calcu- 
lations on the same level of accuracy have been performed 
with the same k-points (in the FBZ), and the volume of the 
"unit cell" has been kept constant. Since MAEs are differ- 
ences in "total energies". Tab. |T] consists of two main blocks, 
namely (left) the evolution of the total energy of the ground 
state, which is the system that is magnetized along the easy 
magnetic axis, here Mll JllOl in agreement with experiment 
(see Refs. [iPllsllellWTllBl^ and theory (see Ref.tH; and 
(right) the evolution in the MAEs with respect to the easy 
axis for the out-of -plane direction, M||[110]. The rightmost 
column shows the converged MAEs for the in-plane direc- 
tion that is perpendicular to the easy axis, i.e., the hard axis 
M||[001]. 

The far left column lists the number of W layers added un- 
derneath the pseudomorphic Fe surface layer. Columns 2-7 
show results for the total energy of the systems with magne- 
tization along the easy axis. Column 2 shows the total en- 
ergy in MeV of the cheapest calculation (lowest value of Tc 
and lowest number of k-points) performed in this study. Col- 
umn 3 shows the difference in the total energy (with respect 
to the results in column 2) when one increases the kinetic en- 
ergy cutoff. One observes that the total energy increases with 
about 340 meV per W atom. In column 4, the density of the 
sampling in k-space is increased, with changes in the total en- 
ergy compared to the results in column 3 of less than a meV 
per atom. Increasing Tc again in column 5, the total energy 
changes at most 100 meV per W atom, and the calculation 
has become so stable that a further increase in the sampling 
in k-space (column 6) gives rise to changes in the total energy 
of less than 1 meV Increasing Tc beyond the value of 9 in- 
troduces ghost bands, and self-consistent minimization of the 
total energy is not possible. With the small changes in the to- 
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tal energy occuring in column 6, however, together with the 
stabiUty seen in the MAEs in the right half of the table, leads 
us to conclude that the calculation has converged with respect 
to Tc and #k. Thus, in column 7 we have listed the total en- 
ergies in meV for the calculation from column 6 to the /jeV- 
level, together with an estimate of the uncertainty obtained in 
the self-consistent minimization procedure. This uncertainty 
is calculated from the total energies of the last three cycles 
of the self-consistent minimization of the total energy. Let n 
refer to the last cycle, then the uncertainty e(«) on the total 
energy V{n) in the last cycle is calculated as: 

e(n) = [|y («) - y(« - 1)1 + \V{n) - V{n - 2)|] /2. (2) 

The smallest uncertainty achieved in the present calculation 
is 2.5 /jRy (34 /veV). We notice in passing that at 7 and 8 
substrate layers, the uncertainties become larger. Beyond 8 
substrate layers, the calculations become less stable, i.e., the 
required level of accuracy in order to calculate well converged 
MAEs could not be reached. 

In the right half of Tab. U are listed the MAEs that come 
out, when the calculations in columns 2-7 are repeated for 
M|| [110] and the results for the easy-axis calculation are sub- 
tracted. Hence, since all "total energies" are negative, positive 
numbers indicate that the absolute value of the total energy is 
smaller than that of the easy-axis calculation, and thus we can 
conclude that the axis is "harder". Similarly, negative values 
indicate that the axis is "easier", and one notices immediately 
that in the two right-most calculations (columns 11 and 12) 
there are no negative values (for 71 = 9). Column 13 lists the 
MAEs for M 1 1 [00 1 ] . Here we have left out the evolution, since 
it shows behaviour similar to the MAEs for M||[110]. It is 
important to underscore the fact that each of the numbers that 
appear in Tab.Uis a result of an individual self-consistent min- 
imization of the total energy and the fluctuations in the charge 
density for the specific configuration. In Tab.|I] three columns 
are in red. The numbers in these columns are the final, con- 
verged, results, and are listed with their individual level of 
uncertainty. Naturally, the uncertainties of the total energy for 
the easy axis are small compared to those of the MAEs. The 
uncertainties listed for the MAEs are calculated as the sums 
of the uncertainties for the total energies of the two calcula- 
tions involved in getting a single MAE. One should bear in 
mind that the uncertainties between the results for M||[110] 
and M 1 1 [001] are not simply an addition of the uncertainties 
listed in columns 12 and 13 of Tab.|I] since the uncertainties 
listed are on the MAEs, not the total energy, and thus with 
respect to M||[liO]. The relative uncertainties between the 
results for M||[110] and M||[001] are (in meV): ±0.20 for 1 
W layer, ±0.07 for 2 W layers, ±0.20 for 3 W layers, ±0.31 
for 4 W layers, ±0.19 for 5 W layers, ±0.16 for 6 W layers, 
±0.39 for 7 W layers, and ±0.37 for 8 W layers. With these 
values in mind, it is safe to conclude that the magnetic hard 
axis is in the [001] -direction. 

In order to visualize the influence of the substrate on the 
magnetocrytalline anisotropy energies, we have in Fig.[T]plot- 
ted the MAEs relative to the easy axis (M||[110]) with re- 
spect to the number of W layers in the substrate under the 
Fe(llO) monolayer, along with their respective uncertainties 
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FIG. 1: (Color online) Magnetic anisotropy energies relative to the 
easy axis (M||[ll0]) are shown as a function of the number of W 
substrate layers underneath the Fe(llO) surface layer. Data are 
taken from Tab. |I] and the text, and the uncertainties have been in- 
cluded as vertical error bars in the plot. The red (full line) and blue 
(short-dashed line) curves show the MAEs of the [001]- and [110]- 
directions, respectively. The green (long-dashed line) curve shows 
the energy difference between samples magnetized along the [001] 
and [110] directions, which also can be called the degree of hardness 
of the hard axis (see text). To guide the eye, cubic splines have been 
used to connect the data points. 



(indicated by a vertical bar at each data point). In order to 
guide the eye, cubic splines have been included to connect 
the data points. The red (upper) line shows the evolution of 
the MAE of the in-plane hard [001] axis, while the blue line 
shows how the MAE of the direction perpendicular to the sur- 
face [110] evolves when adding substrate layers. The green 
(bottom) line shows the MAE between the hard and the per- 
pendicular axes, and can thus be interpreted as the "hardness" 
of the hard axis. As long as it is positive, the hard axis is along 
the in-plane [001] -direction, whereas if it became negative the 
hard axis would shift to being out-of -plane. 

The trend in the values of the MAEs is clear At one layer 
of W, the MAEs are small, but at two layers of W one sees a 
dramatic increase, and the MAEs reaches their maximum val- 
ues. Adding a 3rd layer of W, the MAEs drop again, and from 
the 4th layer and on they get almost constant, in agreement 
with Ref. |25l The value at the "bulk" end (7-8 substrate lay- 
ers) is, with its 3 meV quite close to the experimental value of 
4.2 meV obtained recently by Pratzer, Elmers, et al. (Ref.[35l) 
for Fe monolayer stripes on W(l 10). That the values differ by 
about 1 meV could have its origin in the fact that the exper- 
iment is made on stripes. Since the stripes have lower sym- 
metry than a perfect monolayer we would expect the MAE to 
be slightly higher for the stripes. Older experiment o 1 2, 1 7,36.37 
give values of the MAEs in the range from about 0. 1 1 meV to 
about 6.5 meV, depending on the setup, temperature, and film 
thickness. 

Tableland Fig. [T] reveal a magnetic hard axis along the 
[001] -direction. That the absolute values^ differ from those 
of Ref.[l3|can be attributed to the fact that we are using asym- 



4 



TABLE II: Magnetic anisotropy energies for different magnetization 
directions in 1 ML of Fe on top of 4 ML of W( 110). The anisotropy 
energies are with respect to the value of the total energy for M 1 1 [ 1 1 0] . 
Here, Tc = 9 and #k=2025, and the uncertainties are calculated as in 
Tab.Uby use of Eq. The first four rows are in-plane and the last 
four rows out-of-plane (see text). 



Magnetization direction 


MAE (meV) 


5° from [110] towards [001] 


0.15±0.09 


[111] 


0.83±0.17 


5° from [001] towards [lIO] 


2.82±0.09 


[001] 


2.82±0.19 


5° from [001] towards [110] 


2.53±0.22 


5° from [lIO] towards [110] 


2.31±0.21 


5° from [111] towards [110] 


2.59±0.20 


[110] 


2.34± 0.15 



metric structures (with an Fe layer only on one side of the 
slab) in the present work, whereas Ref. [l3| used symmetric 
structures (with an Fe layer on both sides of the slab). 



B. Magnetic properties of the electronic structure 

In order to verify that the easy axis is really along one of the 
main crystallographic axes (and along the [1 10] -direction) we 
have performed a number of calculations of the total energy 
where the magnetization direction has been shifted by 5° away 
from the main crystallographic axes (in various directions, as 
indicated in Tab.HIli. These calculations have been performed 
for four substrate layers only. Thus, we have in Tab. Ullprinted 
the MAEs with respect to the magnetic easy axis for Tc = 9 
and #k= 2025. In Tab.HIl the first four lines show the devel- 
opment of the MAE as we change the magnetization from the 
easy axis to the hard axis in the surface plane, the next three 
lines show how the MAE changes as we go away from the 
surface, and the last line shows the value perpendicular to the 
surface. The values along the hard and perpendicular axes are 
taken from Tab. U It is important to stress that in order for 
the results in Tab. HI] to be trustworthy (since the MAEs are 
calculated from total energies), one needs to make the calcu- 
lations identical in the sense that only the direction of mag- 
netization is allowed to change. The critical quantity to keep 
identical here is the k-space sampling-mesh (total energies are 
k-dependent, as is evident from Tab.|I]i. The results in Tab.HJ 
are therefore all obtained with the exact same k-points (in the 
Full Brillouin Zone), regardless of the changes in symmetry. 

Of the first four lines in Tab. |II] the first line shows the re- 
sult when the magnetization direction is shifted by 5° in the 
direction of the hard axis. When the magnetization is shifted 
further in the direction of the hard axis, one passes the crys- 
tallographic [lll]-direction (second line of Tab. HHi that could 
be a natural place to look for an easy or hard axis of magneti- 
zation. As we can see, it is neither the easy- nor the hard axis. 
Close to the hard axis (c), the anisotropy at 5° away from it has 
a value that is very close to (equal to, in fact, within our preci- 
sion) that of the hard axis. The value for the hard axis MAE is 
shown for reference in the fourth line. Going 5° out of the sur- 



face from the three axes [110], [111], and [001], we again see 
values between those of the easy and hard axes. In fact, the 
change in the MAE by going 5° out of the surface plane from 
the easy axis is quite strong. Thus, we conclude from Tab.HJ 
that evidence is strong enough for to confirm that the easy axis 
is along the (in-plane) [1 10] -direction, in agreement with pre- 
vious experimenta l'i'^i^'^'^'^'^'i^^ and theoretica l^^i^^ findings. 

Tab. ini] lists the magnetic spin- (jjs) and orbital (jjl) mo- 
ments per atom for the three different magnetization direc- 
tions we consider. First, we notice in general that the mag- 
netic spin moment of the Fe atom is enhanced in comparison 
to the bulk fee Fe value of 2.2 /jB. Second, in all three cases, 
after peaking at about 2.61 /jB for 3 substrate layers, the Fe 
atom takes on a magnetic spin moment of about 2.56 /jB, and 
the first W layer under the Fe surface layer takes on a moment 
of about 0. 1 jjB in the opposite direction (hence, it is antifer- 
romagnetically coupled to the Fe layer), consistent with the 
result of Ref. 39. The remaining magnetic spin moments, in- 
cluding the interstitial moments, are negligible (less than 1.5% 
of the Fe moment). 

Looking at the orbital moments in Tab.lIIIl(^f ). we observe 
that while the spin moments do not change their sign with 
a change of the magnetization axis, the orbital moments do. 
With magnetization along the easy axis (Mj) [1 10]), the orbital 
moment of the Fe surface layer is very small, and in general 
couples antiparallel to its spin moment. For the cases with I- 
2 W substrate layers, the orbital moment in the first substrate 
layer couples antiparallel to its spin moment (and parallel to 
the spin moment of Fe). However, already with the addition 
of the the third substrate layer, the orbital moment of the first 
substrate layer couples in parallel to its spin moment (and an- 
tiparallel to the Fe spin moment). For magnetization along 
the hard axis (M|| [001]) and perpendicular to the surface plane 
(M|| [1 10]), the coupling picture is opposite. The orbital mo- 
ment of the Fe atom is still very small, but now coupled in 
parallel to its spin moment, and the orbital moment of the first 
W substrate layer couples in parallel to its spin moment for the 
cases with 1-2 W substrate layers (and antiparallel to the spin 
moment of the Fe surface layer). The rather large orbital mo- 
ments of W 1 for = 1 and 2 might occur due to the fact that 
these systems tend to have molecular properties rather than 
solid-state ones. They might originate in a combination of 
the large spin-orbit coupling found in W in combination with 
the large Fe moment. From the addition of the third substrate 
layer and onwards, the orbital moments of the first W atom 
couple antiparallel to their spin moments (and parallel to the 
Fe spin moment). Thus, with respect to the orbital moments, 
we conclude (i) that three substrate layers already deliver a 
converged result, and (ii) that in agreement with the observa- 
tions made by Refs J3j,42i Hund's 3rd rule is broken for the 
W substrate. 

For the perpendicular magnetization direction (Mj) [1 10]), 
our results for /j^ differs qualitatively from those obtained in 
Ref. In order to explain this difference, a number of cal- 
culations were performed on the symmetric slab, based on 
the hypotheses that the difference is caused by (i) instabili- 
ties in the computer code, (ii) the perturbative final addition 
of the spin-orbit coupling in Ref. [l3l (iii) the amount of "vac- 
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FIG. 2: (Color online) Evolution in the density of d-states (d-DOS) for an Fe(l 10) monolayer magnetized along the easy axis (M|| [1 lO]) when 
up to 8 W(llO) substrate layers are grown underneath it. The plot to the left (a) shows the density of states for the d-electrons of the Fe atom 
for an Fe monolayer with a single W substrate layer underneath. The d-DOS for the minority spin has been plotted with reverse sign for ease 
of understanding. The plots (b)-(h) shows the changes in the d-DOS of the Fe surface layer when additional W substrate layers are added 
underneath. In order to guide the eye, the scales are kept in constant proportion in all plots. On the energy scales, the Fermi energy is taken as 
the reference point (Ef = 0). In all plots, red (full) lines are the results for the minority spin channel and green (dashed) lines are the results 
for the majority spin channel. In order to assure a good view of the plots we have moved the scales to the border of the figure. 



TABLE III: Magnetic spin (jjg) and orbital (jji) moments per atom (in Bohr magnetons, /jb) are shown for the pseudomorphic Fe surface layer 
as well as the first W substrate layer (Wi) as a function of the number of substrate layers (AV) when the magnetization is along the in-plane 
easy axis (columns 2-5), the in-plane hard axis (columns 6-9), and perpendicular to the surface (columns 10-13). 

Easy axis : M 1 1 [ 1 1 0] Hard axis : M 1 1 [00 1 ] Perpendicular: M 1 1 [ 1 1 0] 

Surface: Fe Substrate: Wi Surface: Fe Substrate: Wj Surface: Fe Substrate: Wi 



Nw 




ML 


Ms 


ML 


Ms 


ML 


Ms 


ML 


Ms 


ML 


Ms 


ML 


1 


2.562 


-0.00023 


-0.135 


5.792 


2.561 


0.00005 


-0.136 


-2.184 


2.560 


0.00007 


-0.135 


-0.584 


2 


2.574 


-0.00010 


-0.100 


2.517 


2.574 


0.00009 


-0.099 


-3.449 


2.570 


0.00007 


-0.102 


-1.359 


3 


2.616 


-0.00006 


-0.089 


-0.021 


2.613 


0.00005 


-0.091 


0.021 


2.609 


0.00005 


-0.094 


0.024 


4 


2.561 


-0.00006 


-0.104 


-0.025 


2.557 


0.00005 


-0.105 


0.024 


2.559 


0.00005 


-0.104 


0.028 


5 


2.567 


-0.00006 


-0.098 


-0.024 


2.563 


0.00004 


-0.100 


0.023 


2.565 


0.00005 


-0.099 


0.026 


6 


2.568 


-0.00005 


-0.101 


-0.024 


2.567 


0.00004 


-0.101 


0.023 


2.567 


0.00005 


-0.100 


0.027 


7 


2.559 


-0.00005 


-0.104 


-0.025 


2.558 


0.00004 


-0.104 


0.024 


2.558 


0.00005 


-0.104 


0.028 


8 


2.565 


-0.00005 


-0.103 


-0.025 


2.563 


0.00004 


-0.103 


0.024 


2.562 


0.00005 


-0.103 


0.028 



uum" between the slabs in the supercell calculation, or (iv) 
the breaking of the symmetry in the slab (removing one of the 
Fe surfaces). A repetition of the calculation in Ref. [Tsl elimi- 
nates the first hypothesis. A comparison between the repeated 
calculation and one with spin-orbit coupling included selfcon- 
sistently leads to elimination of the second hypothesis.'*" In 
the symmetric slab, the Fe layers may couple electronically 
to each other either through the W substrate or through the 
vacuum. In Ref. the distance through the W substrate 
is 24.52 Bohr, and the distance through the vacuum is 26.01 
Bohr. Thus, in order to test hypothesis (iii) against hypothe- 
sis (iv) we increase the distance between the two Fe layers on 
the vacuum side to 76.54 Bohr^' The result of the calculation 
with the increased spacing on the vacuum side between the 
two Fe surfaces is [in support of hypothesis (iii)] an orbital 
moment for the Fe layer of 0.00005/jb and an orbital moment 
for the first W layer of 0.028a(b, in agreement with Tab. [Till 
Thus, in agreement with Refs I42ll43i we may conclude from 



this little excercise on the symmetric slabs and from Tab. |III] 
that both the size as well as the alignment of the orbital mo- 
ments in the Fe surface layer and the first W substrate layer 
are not only very sensitive to the local structure, but also to 
the direction of magnetization in the Fe layer 

In order to determine the influence of the substrate layers 
on the electronic properties of the magnetic surface, the evo- 
lution in the density of states for the d-electrons (d-DOS) of 
the Fe atom has been plotted in Fig. |2] (easy-axis magnetiza- 
tion, M||[liO]), beginning with the Fe/W(110) bilayer resuh 
in Fig. 12 a). In Figs.|2jb)-(h), the differences due to the addi- 
tion of a second, third, fourth, fifth, sixth, seventh, and eighth 
substrate layer, respectively, are depicted on the same scale as 
the d-DOS of the bilayer. It is clear from Fig.|2]that the main 
features of the d-DOS of the Fe(llO) surface layer are to be 
found in a band between 4 eV above and 5 e V below the Fermi 
energy {Ef = in Fig.O. Also, we observe from Fig.|2]that 
by adding substrate layers, the changes to the d-DOS of the Fe 
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FIG. 3: (Color online) The density of d-electrons of the Fe atom (d- 
DOS) is shown for the calculation that has converged with respect to 
the number of substrate layers, i.e., for Nv/ = 4. Results are shown in 
the combined plot (a) for the easy (red, full line), hard [green, long- 
dashed line), and perpendicular (blue, short-dashed line) directions 
of magnetization, and (b)-(c) show the differences in the Fe d-DOS 
between (b) the easy and hard axes and between (c) the easy and per- 
pendicular axes, respectively. In (b) and (c), the differences in the 
minority spin channel are plotted in red (full lines), and the differ- 
ences in the majority spin channel in green (dashed lines). For ease 
of understanding, the d-DOS of the minority spin has been plotted 
with reverse sign, as in Fig. [2] Again, the zero on the energy scale 
corresponds to the Fermi energy and, again, we have put the scales 
to the border of the figure. 

surface layer become less and less pronounced, as the number 
of substrate layers go up. Since already after the fourth or 
fifth substrate layer, the d-DOS of the Fe surface layer has 
converged to within a few percent of what it would be on bulk 
W(llO), the d-DOS supports the conclusion that the MAEs 
have reached their bulk value after adding 4-5 substrate lay- 
ers. 

When the d-DOS has converged, the differences in d-DOS 
due to changes in the magnetization direction are very sub- 
tle. In order to illustrate this, we have in Fig. [3ja) plotted 
the d-DOS of an Fe(llO) monolayer with four W substrate 
layers underneath it for the three cases where the magnetiza- 
tion direction is along one of the three main crystallographic 
axes. From Fig. ^a.) we notice that the d-DOS is now al- 
most identical in all three cases. In order to explore the sub- 



tle differences in the d-DOS between the different magneti- 
zation directions, we have in Figs. [3jb)-(c) additionally plot- 
ted the differences in the d-DOS between the easy axis, the 
hard axis, and the axis perpendicular to the surface, i.e., the 
two quantities d-DOS(M||[110])-d-DOS(M||[001]), and d- 
DOS(M||[ll0])-d-DOS(M||[110]), respectively. In contrast 
to the evolution of the d-DOS when adding substrate layers, 
the differences in the d-DOS between the different magnetiza- 
tion directions are too small to be plotted on the same scale as 
the d-DOS itself. Since the MAEs are in the meV-range, this 
is to be expected, since large changes in the d-DOS would 
lead to large MAEs. 

IV. CONCLUSIONS 

We have shown that the magnetocrystalline anisotropy en- 
ergy in an Fe(llO) monolayer on W(llO) can be converged 
with respect to the thickness of the substrate using ab ini- 
tio methods. After showing large changes for the first few 
substrate layers, it stabilized close to 3 meV, close to the ex- 
perimental value. In addition, the directions of the easy and 
hard axes came out consistently, and in-plane. As expected 
from the large variations in the MAEs with respect to the ad- 
dition of the first few substrate layers, also the density of the 
Fe d-states vary a lot during the addition of the first few sub- 
strate layers, after which it stabilizes. At the bulk-like sub- 
strate thicknesses, the differences between the density of d- 
states between the different magnetization directions reflects 
the fact that the MAE is in the meV regime (they are very 
small). 
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